Decay constants in the heavy quark limit 
in models a la Bakamjian and Thomas 



V. MorenasQ 

Laboratoire de Physique Corpusculaire 
Universit Blaise Pascal - CNRS/IN2P3 F-63000 Aubiere Cedex, France 

A. Le YaouancQ, L. Oliveif], O. PeneQ and J.-C. Raynal 

Laboratoire de Physique Theorique et Hautes Energies^ 
Universite de Paris XI, Batiment 210, F-91405 Orsay Cedex, France 



Abstract 

In quark models a la Bakamjian and Thomas, that yield covariance and 
Isgur-Wise scaling of form factors in the heavy quark limit, we compute the 
decay constants and of S-wave and P-wave mesons composed of 

heavy and light quarks. Heavy quark limit scaling \f~M f = Cst. is obtained, 
and it is shown that this class of models satisfies the sum rules involving 
decay constants and Isgur-Wise functions recently formulated by us in the 
heavy quark limit of QCD. Moreover, the model also satisfies the selection 

(n) 

rules of the type = that must hold in this limit. We discuss different 
Ansatze for the dynamics of the mass operator at rest. For non-relativistic 

kinetic energies - — the decay constants are finite even if the potential V(r) 

has a Coulomb part. For the relativistic form \J-p 2 + m 2 , the S-wave decay 
constants diverge if there is a Coulomb singularity. Using phenomenological 
models of the spectrum with relativistic kinetic energy and regularized short 
distance part (Godfrey-Isgur model or Richardson potential of Colangelo et 
al.), that yield p 2 = 1 for the elastic Isgur-Wise function, we compute the 
decay constants in the heavy quark limit, and obtain fg — 300 MeV, of the 
same order although slightly smaller than in the static limit of lattice QCD. 

We find the decay constants of D** with j = — of the same order of magnitude. 
The convergence of the heavy quark limit sum rules is also studied. 
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1 Introduction 



In a recent paper jl| we have shown that quark models of hadrons with a 

fixed number of constituents, based on the Bakamjian-Thomas (BT) formalism ||, 

yield form factors that are covariant and satisfy Isgur-Wise (IW) scaling || in the 

heavy mass limit for one of the quarks. In this class of models a lower bound is 

3 

predicted for the slope of the heavy meson elastic IW function p 2 > - . Moreover, 
the model satisfies the Bjorken-Isgur-Wise sum rule 0] that relates the slope of the 
IW function to the P-wave form factors Ti/2(w), t 3 / 2 (w) at zero recoil 0. We have 
also explicitly computed the P-wave meson wave functions and the corresponding 
inelastic IW functions || within the model, we have made a numerical study of p 2 
in a wide class of models of the meson spectrum (each of them characterized by 
an Ansatz for the mass operator M, i.e. the dynamics of the system at rest) 0, 
and a phenomenological study || of the elastic and inelastic IW functions and the 

corresponding — for B — ► D, D*, D**£v. 
aw 

The purpose of the present paper is to study the decay constants of heavy mesons 
within the same approach. 

2 Decay constants in the B-T scheme : heavy 
quark scaling 

We want to compute in the model the decay constants of mesons Qq, where Q 
and q are a heavy and a light quark, defined by 

< P(0-)(v)\Af\0 >=N f P ^M Vfl 

< V{V)(v,e)\Vf\Q >=Nf v VM e J (1) 



for S states, where iV = 




and M is the mass of the corresponding bound state 



(M = M P or M = M v ), and similarly : 



< 1/2 O+^I^IO >= N fW*>VM v» 

< 1/2 \+(v,e)\Af\0 >= N gWy/M e£ 

< 3/2 l+(v,e)\Af\0 >= N gWy/Me^ (2) 
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for P states. The model satisfies the heavy quark limit relations (we use the phase 
convention for states of ref. 0) 



fp = fv = f 



(3) 



g d/2) = _ /( l/2) 0(3/2) =Q (4) 

In terms of the internal wave function we get (the superindex (n) means any radial 
excitation), for 5- wave mesons : 



Mf (n) = JN C V2 



dp 2 1 



V '(P2 ■ v)(p 2 ■ v + m 2 ) ip {n) (k 2 ) 



(2tt)3 p o 

where the wave functions at rest of 0~ mesons are given by 



(5) 



and for P wave mesons 



dp 



^(p)| 2 = l 



(27T) 



(6) 



dp2 1 y/P2 ■ V 



^i/ 2 (k 2 ) |(P2 -v) 2 -m 



(2tt) 3 ^ y/p 2 -v + m 2 
where the wave function at rest of + mesons is given by 



(7) 



In these expressions, 



i 

dp p 2 
(27p 



cr ■ k 2 ) cr' 



S1,S2 



yl^i/ 2 (P 



(8) 



p 2 = B v k 2 



(9) 



where _B„ is the boost operator. The expressions (|5|) and ([/]) are Lorentz invariant 
since, due to the rotational invariance of y?( n )(k 2 ), (pjj 2 (k. 2 ), they are functions of 



1/2 

f 2 = 1. One can take simply the rest frame v = (1,0). 

Let us give a brief outline of the calculation. Within the model [||, the internal 
wave function in motion of a S-w&ve meson with center-of-mass momentum P (n 
denotes the radial excitation) is given in terms of : 
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*s 2 (p-p 2 , P2 )= x ^ n t=e n [A(R,.)]^^(k 



(10) 



where (k 2 ) is the internal wave function at rest ([5]), and Rj are the Wigner 



SlS 2 

rotations : 



Rj - Bp/ B s P , Bfc. . 
The qQ to vacuum matrix element of a current, e.g. J = 7^75, reads 



(11) 



S1S2 (P - p 2 ,P2) =< 0|J|P -p 2 si;p 2 s 2 >= [m Si (P - p 2 )]*i7 2 7 7 M 7 5 M S2 (p 2 ) (12) 
The matrix element of interest to us 



< 0\A»\P^ >= ^ / 7^ E <W P - P2> P 2 )^S 1 (P - P2, Pa) (13) 

J {Z7T) {*} 



can be written, after some algebra, using (]ToD - (]T2|) : 



1+7° 1 1+7° 
— — — B"S M B 

Pi ' P2 



2 ^0 V PI V ^2 

BlrBj^ B fc2 Bj^B^ B P1 



mi m 2 fnVl . 



P2 



1+7° 



(14) 



and since the projector commutes with the Wigner rotations (PD, one can 



rewrite : 



where u 



< 0|A"|PW >=\fN c f ^ TT 



^P 2 M TT & & /^)(k 2 ) 
1 m I l + V B l + f D , n , 



^ 7"B„B fe 2 



2 \/Pi 

Br/B 



.0 v pi v 



fcl u 



Pi +P2 
(Pi +P2Y 



-.. Using the identities : 



R R 1+ ^ R~1 
t5 u t}fc 2 ^ B u 



m 2 + p 2 l + jx 



2m 2 (k° + m 2 ) 2 



B, 



l + 7 -o-i-o-i _ 1 +^ mi+A 



2 J2mi (&° + mi) 



(15) 



(16) 
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one obtains 



<0 |^>^/^ (P2 , P) 

Tr [ 7 " (m 2 + i> 2 ) (1 + ii) (mj + &)] ^»>(k 2 ) 



with k 2 = B n x p and 



F(p 2 ,P) = ^2^— ^ 



Finally, in the heavy mass limit 



m 2 



(17) 



;is) 



u — > v 
mi 

one obtains : 



— -> w 
mi 



(19) 



< 0|^|P (n) >= 



dp 2 1 ^Vi ■ v 



^2^2 J {2n) 3 p° 2 ^p 2 ■ v + m 2 
Tr[ 7 ^(m 2 + ^)(l + ^)]^")(k 2 ) (20) 

where k 2 = B~ x p2 and <^ n )(k 2 ) depends on k| = (p 2 -w) 2 — m 2 because of rotational 



mvanance. 



Thus, the integrand must be proportional to t> M , yielding equation (|5|) for 
the decay constant. This expression exhibits the expected heavy quark limit scaling. 
Moreover, varying the current, one can easily show relations @. It is also worth 
noticing that expression (f|) becomes the well-known non-relativistic expression || 
if we assume that the light quark is also non-relativistic. 

Let us give some details of the calculation of (§) and the corresponding relations 
(H). Mutatis mutandis , we consider the matrix element < 0|V M |0 + > with the + 
wave functions given by ([8]). After some algebra, the matrix element equivalent to 
(|15D is in the present case : 



< 0IV10+ >= 



dp 2 



11 rR\l p o y p° ^va) 



k 2 ) 



B; 1 B u ( CT .k 2 )B; 1 B u i±^B fc -X 1 ) • (21) 
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Using ([U]) and taking the limit fll~9l) , we obtain 



< o\w >= ^54= f -^1 y f)(fa) 



Tr [-y (ma + (1 + # (-# 2 + P2 • u) (1 - 



(22) 



and using the fact that this expression must be proportional to v^, we obtain expres- 
sion fl7|), scaling invariant in the heavy quark limit. Moreover, varying the current 
and considering instead the matrix element < 0[v4.' z | ljy 2 > one can easily verify the 
first relation (f|). 

The result g( 3 / 2 ) = 0, that must hold in the heavy quark limit deserves a few 
details. The matrix element to be considered is (m denotes the polarization of the 
state 1 + , and n the radial excitation) : 



< 0|^|l+ 2 ,m >= \/N ( 



dp 2 



Ep^ „ Jkf j mi jm 2 



(2tt)^M ^ 



P° V Pi V ^2 



Tr { ( ^a 2 B- i 1 z 7 ° 757 ° 7 - 75 B P2 i±^') [D(R 2 ) 



^ (n) (k 2 ,m) 



^(Ri) 



(23) 



where the wave function is, for the l + j = — states | 



s (m) 



-k 2 ^= (k 2 x cr) 

3 VG 



(T2 <4/ 2 (k 2 ) 



(24) 



where is a unit vector, and the rotational invariant function <^3™ 2 (k 2 ) is norma- 
lized according to (|S|). The equivalent expression to (^Tj) will be 



dp 



Tr 7 "B U B 



i + 7 ° 



.M 



^k 2 + -j=(k 2 x<r) 



We need to compute the expression B„e^ • y |k 2 + ^g(k 2 x cr) B^ 1 in the heavy 
quark limit. After some algebra, one gets : 



B„ e (m) 



T^ 2 x a) 



b: 1 - -4= 



i 



4^ • P2 



(26) 



6 



where 4 m) is the axial meson polarization in the heavy quark limit. Using ((T6|) in 
the limit (|19D we obtain finally, 



W c 1 



' ' 3/2 8^^^ (2n) 3 p% y/p 2 ■ v + m 2 * 3/2y V 

Tr (ma + i> 2 ) (1 + [~t ] i>2 + ^V(P2 • «) - 4 m) " ft] (1 + 0} (27) 

This expression is covariant and satisfies heavy quark scaling. In the rest frame, 
particularizing to m — 0, one gets, because of rotational invariance of ip 3 / 2 (k 2 ) : 



^ 1 



< 0U*|1+ ,P = 0,m = >- _ , , n 



dk 2 1 



= ^3/2 (k 2 



2 + m 2 



{-[(^) 2 + (^) 2 ]+ 2 ^) 2 } = • (28) 

Therefore the second relation (|J), g( 3 / 2 ) = 0, follows. The vanishing of g( 3 / 2 ) 
can be seen also from the following covariant argument. Contracting (|27|) with the 
four-vector 8)™' we see that 

0(3/2) is 

proportional to the integral 
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(3/2) 



C/P 2 1 y/p 2 ■ V 



y/p 2 ■v + m 2 
Two types of integrals appear : 



V?3/ 2 (k 2 ) Wr ) -P2? + {v-p2) 2 -m 



(29) 



L flU 



dp2 1 y/P2 ■ V 

(27r) 3 ^ y/p 2 ■ v + m 2 

dp 2 1 y/p 2 ■ V 



<^3/2(k 2 ) P2p, V2v 



(2-kYpI y/p 2 -v + m 2 
From covariance it follows that 



^3/2 (k 2 



(30) 
(31) 



I = C 

where A, B, C are constants. Contracting with g^ v it follows 



(32) 



45 + A = m\C 



and therefore 



(33) 



^ ~ AB + A - mlC = 



(34) 
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3 Sum rules 



Let us now show that the QCD heavy quark limit sum rules [[K]] 

n J 

Tx/2(w) ^E^n ( /lH = i (36) 

are satisfied within the present B-T scheme. 

In order to prove Bjorken sum rule, in ref. || we used the completeness relation 
at fixed P that holds in the B-T formalism for the internal wave functions. We 
proceed here in the same way. We need to compute expressions of the form : 

E<0|O|P',nxP>|O|P,0> (37) 

n 

where < 0|O|P', n > and < P', n|0|P, > will be related respectively to the decay 
constants and to the Isgur-Wise functions. We need P' different from P in order to 
demonstrate the sum rules (|35l)-(j36l) for any value of the scaling variable w. 

To obtain the sum rule ( p5|) for S-waves, the simplest way is to take P' = for 
the intermediate states and choose the currents O = 7°7 I/ , and O = 7°7 At 75 with 
/i = 0. Then, only the 0~ intermediate states contribute, because the 1 + states 
do not at P' = 0, since then the time component of the polarization Eq = 0. 
Alternatively, to obtain the sum rule (|36|) for P-waves, we will take P' = and 
O = 7°7 y 75, and O = 7°7 At with fi = 0. Then, only the + intermediate states 
contribute, since the 1~ do not at P' = 0. 

We start from the completeness relation at fixed P' : 

E (P' " P2, P2) *S ( P ' - ^ P^) = 5 « ( 27r ) 3 ^ - P2) • (38) 

n 

We need the matrix elements 



< P',n|O|P,0 >= 




CS ) *(Pi>P2)a; Sl (p;,Pi)v&S(pi,p 2 )(2vr) 3 5(p 2 - p' 2 )^ S2 (39) 
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< 0|O|P,n >= JN t 



dp 2 

(2tt) 3 



^0 SlS2 ( Pl ,p 2 )^( Pl ,p 2 ) 

K} 



(40) 



related respectively to the IW functions and decay constants. In these expressions 



0(p'i,Pi) 



VmiT^i +7 o 1 + 7° 

2 p 'i Pl 2 



PiPx 



0(P -p 2 ,p 2 ) = W^o"^— ^B p ,« 7 7 5 OB P2 ^— 



(41) 



(42) 



and here O and O are the Dirac matrices of the corresponding currents. After some 
algebra we obtain, using the completeness relation (38) : 



< 0|O|P',n >< P',ti|O|P,0 >= 



(t^EE S[S2 (P' - pa, p 2 ) O s , Sl (P' - pa, P - p 2 ) ^ (P - p 2) p 2 ) 



(43) 



where the ground state wave function is given by [Fl j 



n 



D(R 1 )^(°)(k 2 ) J D i (R 2 ) 



N M ° i=l,2 

and using @, in an obvious notation : 



(44) 



Y, < 0|O|P',n >< P',n|O|P,0 > 



C?P2 T-r M n i l (0)/, \ 



Tr [o* (P' - p 2 , p 2 ) O (P' - p 2 , P - p 2 ) D(R 1 )o- 2J D*(R 2 

where the matrices R; are given by ([TT|). 

To isolate the 0~ states, we take, as argued above : 



(45) 
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< 0\A°\O,n >< O,n\V"\P,0> 



dp 2 



(2tt)3\ Mo U ^ V rfpf V V P° 



Tr 



1 + 7° 



2 a 2 B-z 7 u B 



'Pi 

1 + 7° 

P 2 2 



i±^B-/ T °rB Pl i±^ 
2 Px 2 



(46) 

D(R!)(t 2j D'(R 2 ) 



Using a 2 (B p )V 2 = B" 1 , a 2 {^fa 2 = 7 ", B p } ^ , ? , . ^ , . 



H , the fact 



that -(1 + 7 ) commutes with rotations, and the relations 
2 



B„B fcl B v 1 



mi + 



^2mi(fcJ + mi) 



m 2 + fbfo 



: (47) 

'2m 2 (fc§ + m 2 ) 

one obtains, after some algebra, in the heavy quark limit (pi — > mit>, — > mi), for 
the r.h.s. of (ESI 



" L 1 /" rf P 2 1 / — /To 1 (0)/, \ 

Tr (i + 7°) 7^(1 + flB^B^B; 1 } = 

- 1 i y^J _ l (0)( } 1 

V V4vV°\/2Vo^ (2tt) 3 ^ ^■t; + m 2 y/T v 7 2 

TY {(1 + y) Yd + « (m 2 + «} = ^^VM/'") + V) . 



On the other hand, the l.h.s. of (HGI) reads : 



< 0|A°|0,n >< 0,n\V v \P,0> 

n 

l : 



^^^(^(f + /) (49) 



and therefore the sum rule ( p5|) follows. 

The sum rule for the P-states (^6|) follows straightforwardly in a similar manner 
by taking P' = and O = 7°7 iy 75, O = 7°7 At with /i = 0, since then only the 
+ intermediate states contribute, as pointed out above. However, to illustrate the 
methods of calculation in the BT formalism, we will now verify the further sum rule 
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( p6|) from the general expressions for f^l and t^Mw), without appealing to any 
particular frame. 

We want to evaluate the expression 

Ti /2 H/ (0) = EA ( /1^H • (50) 

n 

From the explicit expressions (|7|) for fy\ and Ty\ (w) from ref . || : 
1/21 2^3-/ (2ir) 3 p? s f(p^. v + m2 )(pf 2 . V ' + ma) V 



(p' 2 ■ f / )(p' 2 .f + m 2 ) — (p 2 " V )(P2 " u + wm 2 ) + (1 — w)m 2 



1 — u> 

with 



(51) 



fc£ = B-y 2 ^ = b;,v 2 . (52) 

The computation of the sum (|50|) leads to the expression 

E?$04)V$(P2) = 67r 2 ^5(|p 2 | - |k 2 |) (53) 

n P2 K 2 

where the r.h.s. follows from ref. ||. One can then perform the integration over p 2 
that amounts to replace all p 2 by (p 2 ■ v). Realizing then that 



(p 2 ■ v')(p' 2 ■ v + m 2 ) — (p' 2 ■ v)(p 2 ■ v + wm 2 ) + (1 — w)m 2 

1 — w 

k / 2 (_ x + P' 2 ■ v' - w(p> 2 ■ v) \ 
2 V (Pz ■ v ~ m a)(l -w)J 



one gets 



2T 1/2 H/(°)=2E/^H1' 



It; 



/( o) Wc r d P ' 2 1 yc^v) ^^ pj^jjj _ (55) 



7T J (27r)3p 2 ° y/fa-v 1 + 702) ' ' 1_u ' 
The second term involves the integral 

rfp' 2 1 V (P2 • «0 {0) , , 



(27r) 3 p 2 ^(p' 2 . t; / + m2 ) 
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where C is some constant, because of covariance. Inserting the r.h.s. of ([56]) into 
( [55D we see that the second term of the r.h.s. of ( p5|) vanishes for w ^ 1, and then 
the sum rule ( |36"D follows. For w = 1, one can choose v = (t>°,v), v' = (1,0) and 
make an expansion of (w — 1) _1 for v — > 0. Then, the second term in the r.h.s. of 
( p5| ) vanishes from rotational invariance. The same method allows also to obtain 
the sum rule fl35|) along similar lines. 

4 Numerical results 

As explained at length in ref. [I], the dynamics in the B-T formalism depends 
on the form of the mass operator M at P = 0. This mass operator can be of any 
form 

M = K({k l })+V({r l , Pl }) (57) 

and one obtains covariance of the form factors and IW scaling with only the very gen- 
eral assumption of rotational invariance of M. We are going now to give the results 
for the decay constants for various Ansatze of the operator M, not only various 

forms for the potential V^({rj,pj}), but also for the kinetic energy K({ki}), that 

k 2 / 

can be taken to be of the non-relativistic — — or relativistic wk 2 + mf forms. We 

2m,i v 

choose such models in order to emphasize the physics involved in the decay constants, 
sensitive to the short distance part of the potential and to the scaling behaviour of 
the kinetic energy (quadratic or linear in k). Of course, any scheme that we adopt 
should give a reasonable fit of the whole meson spectrum. The success of quark 
models in the description of the spectrum with either non-relativistic or relativistic 
kinetic energies shows that the spectrum by itself does not constrain the form of this 
kinetic energy. Interestingly, we have shown in ref. ]7) that quark models of form 
factors in the BT formalism show a clear preference for the relativistic form of the 
kinetic energy, since they give a slope of the elastic IW function p 2 = 1, while models 
with a non-relativistic kinetic energy give a larger, phenomenologically unacceptable 
value. 

We list here a number of phenomenological quark models of the hadron spectrum 
that we will use, specifying their interesting features : 

1) Isgur, Scora, Grinstein and Wise (ISGW) spectroscopic model [|TT|j (to be dis- 
tinguished from the ISGW non-relativistic model of form factors ) : non-relativistic 
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k 2 

kinetic energy — — with linear plus Coulomb potential. 
2m,- 



2) Veseli and Dunietz (VD) model |T2| : relativistic kinetic energy -y/k 2 + m 



with linear plus Coulomb potential. The Coulomb part is not regularized. 



3) Godfrey and Isgur model (GI) [[L3|] : relativistic kinetic energy y k? + mf with 
linear plus a regularized short distance part. This scheme also incorporates the fine 
structure of the potential. 

4) Richardson potential with relativistic kinetic energy, as used by Colangelo, 
Nardulli and Pietroni (CNP) |14|]. This model exhibits the asymptotic freedom be- 
haviour at short distances, but the Coulomb singularity (logarithmically corrected) 
is regularized by a cut at small r. 

These models are solved numerically using a harmonic oscillator basis. In Table 1 
we give the masses and decay constants of the states n = 0, I = or 1 as a function 
of N max in the different models. N max means the number of radial excitations 
included in the truncated basis, the ground state gaussian plus up to N max radially 
excitated harmonic oscillator wave functions : N max + 1 is then the dimension of 
the truncated Hilbert space. In Table 2 we give the decay constants of a number of 
radial excitations in the GI model. 

Let us begin our discussion with the models with relativistic kinetic energy. A 
first important remark to make is that dimensional analysis implies that a Hamilto- 
nian (or mass operator in the B-T formalism) with kinetic energy of the relativistic 



form y k 2 + mf and a Coulomb potential implies divergent wave functions at the 
origin \I/( n )(0) (or decay constants f^ n ') for 5-waves. The reason for this behaviour 
is that the kinetic and the potential energies exhibit the same scaling properties at 
large momentum or small distances. The relativistic kinetic energy is not efficient 
enough in smoothing the r-space wave function. This divergence of the wave func- 
tion at the origin is the cause of the existence of a critical coupling a crtt in this 
class of models : for a > a CTtt arises the so-called phenomenon of fall in the center. 



For a discussion, see for example the paper by Hardekopf and Sucher ||15|| . On the 
contrary, the decay constants for P-wave mesons remain finite even in this case. 
A short distance Coulomb part corrected by asymptotic freedom exhibits the same 
phenomenon, although /w diverges only logarithmically in this latter case, instead 
of as a power in the former. These general features are exhibited by the model of 



Veseli and Dunietz [12], as shown in the Table, and have been underlined by these 



authors. In our Table, the finite values obtained for a/M /(°) using the singular 
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model for N max = 10, 15 or 20 are just an artifact of the truncation method. 



The GI model jy| , a model of the meson spectrum for all qq, Qq and QQ systems, 
chooses a short distance part with a fully regularized Coulomb singularity. Within 
this model it is possible to compute both types of decay constants and /]! In 
the Table we give the scaling invariant quantities y/Mf^ and \/Mffj 2 for S and 
P-wave n = mesons, computed in the heavy quark limit. 

We have also made the calculation for the CNP model (Richardson potential 
plus relativistic kinetic energy) fll4"f . 

k 2 

In the case of a kinetic energy of the non-relativistic form the S"-wave decay 
constants are finite even in the presence of a Coulomb singularity. This feature is 
examplified by the ISGW |Ll| model in Table 1. Moreover, the decay constants 
are smaller in this case than for the models of relativistic kinetic energy, another 
manifestation of the singular behavior of the latter. 

It is interesting to notice that in the heavy quark limit and for N max = 20 
one gets, for the models with relativistic kinetic energy and regularized Coulomb 
singularity : 



\fMf<® = 0.67 GeV 3/2 (GI model) ; 0.83 GeV 3/2 (CNP model) 
y/Mfl°\ = 0.64 GeV 3/2 (GI model) ; 0.70 GeV 3/2 (CNP model) . (58) 

Applying this asymptotic result to the B meson, one obtains 

f B = 300 MeV - 350 MeV . (59) 

This value is not far away, although slightly larger than the values obtained in 
lattice QCD in the static limit , to which it should naturally be compared, which 



ranges between 220 and 290 MeV [JIB] . However, one should keep in mind that the 
lattice QCD result includes logarithmic corrections absent in our phenomenological 
scheme. 

It is worth noticing that we obtain the same order of magnitude for the decay 
constant fy 2 . This is phenomenologically important because, reasoning within the 
factorization assumption, this means that the emission of D**(0 + ) and D**(l + , j = 

-) is expected to be important in B decays. On the contrary, the emission of 
3 

D**(j = — ) will be suppressed. Table 2 shows that the decay constants of radial 
excitations are of the same order of magnitude as in the ground state. However, the 
error due to the truncation is larger as n increases. 
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Finally, let us study the convergence of the sum rules (|35|) and (|36|) , that we have 
shown formally to hold, in models that give finite results (GI, CNP and ISGW) for 
which we can compute the decay constants f^ n \ fy\ and the IW functions ^ n \w), 
Tyl(w) ||. The convergence of the Bjorken-Isgur-Wise sum rule || has been studied 
in ref. ||. 

Let us define 

N fin) N f( n ) 

* (A °M = E7# (n) H T i/?H = E^| T 5H • (so) 

n=0 J n=0 J 

We compute the sums for N = N max in the different models for various values of 
w, and see how they compare to the r.h.s. of the sum rules ( |3~5"D and fl36l). Let 
us recall that N max is the maximal number of radial excitations included in the 
truncation method (the dimension of the variational base is N max + 1). We show 
the results for the Godfrey-Isgur model in Figures 1 and 2. The Ox axis represents 
^-/{Nmax + 1)- We observe that these sums converge fairly well towards the r.h.s. of 
the sum rules (respectively 1 and \) as we increase N max . For fixed N max we can 

ask how the partial sums IW(id) and T$(w) (N < N max ) behave clS 8b function 
of N, i.e. how fast X^ N \w), T$(w) approach X( Nma "\w), T^ 2 max \w) when N 
increases (N > 0). Let us give the example N max = 20. The convergence is rather 
fast, but it degrades as w increases. Concerning the X sum rule, for w = 1 one has 
trivially X^(l) = X( N ™*\l), because the ground state saturates the sum, since 
£( n )(l) = (n > 0). As one increases w, one one needs to sum up to iV = 3, 4, 5, 6 
or 7 respectively for w = 1.1, 1.2, 1.3, 1.4 or 1.5 to approach X^ Nmax \w) at the 5 % 
level. Concerning the T1/2 sum rule, one needs iV = 3, 4, 5, 6, 7 or 8 respectively 
for w = 1.1, 1.2, 1.3, 1.4 or 1.5 to approach T^ 2 max \w) at the 10 % level. 

For the CNP model and the non-relativistic ISGW model the convergences both 
in N max and in iV for fixed N max are not as good. For the case of the VD model 
one gets finite fy\ and Tyl(w), but the convergence toward the r.h.s. of the sum 
rule (|30|) does not improve as one increases N max . This results from the divergence 
of the denominator /(°) in (|35|), (|36|) when N max — > 00. 

In conclusion, we have studied the decay constants of heavy-light mesons in the 
heavy mass limit in a class of models a la Bakamjian and Thomas, with different 
Ansatze for the dynamics at rest, with non-relativistic or relativistic kinetic ener- 
gies. Each particular model gives an acceptable phenomenological description of 
the spectrum. The models with relativistic kinetic energy, that yield a slope of the 
elastic IW function p 2 = 1 (as shown in ref. J?]), give finite decay constants if the 
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Coulomb singularity of the potential is regularized, as in the GI model. At the B 
mass, one finds Jb slightly larger than in the static limit of lattice QCD. The decay 
constants of D** with j : — - are of the same order of magnitude. Moreover, we have 
shown that heavy quark limit sum rules involving decay constants [10] are satisfied 
by these class of quark models a la Bakamjian and Thomas, and in the case of the 
Godfrey-Isgur model the convergence of the sum rules is quite fast. 
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Table captions 

Table 1. Decay constants of n = 0, I = and n — 0, I — 1, j — — mesons in the 
various models. The ISGW is non-relativistic. In the VD, GI and CNP models, the 
kinetic energy is relativistic. In the VD model the S'-wave decay constants diverge 
due to the Coulomb singularity In the GI and CNP models the Coulomb singularity 
is regularized and the decay constants are finite. N max stands for the number of 
radial excitations included in the truncated variational harmonic oscillator basis. 

Table 2. Decay constants and fy\ for the first radial excitations of £ = and 

£ — 1, j — — mesons. The error in parenthesis is estimated by comparing the number 
of radial excitations included in the truncated variational harmonic oscillator basis 
N m ax = 20 and N max = 10. For n > 5 the error becomes larger than 20 %. 



Figure captions 



Figure 1. Convergence of the heavy quark limit S-wave sum rule (|35|) in the GI 
model [O] for different values of the scaling variable w as one increases N max , the 



number of radial levels included in the truncated variational harmonic oscillator 

N fi n ) 

basis. The Ox axis represents (N max + 1) and the Oy-axis "ttfvTs n { w )> that 

n=o p u > 

the sum rule predicts to be equal to 1 for N max — > oo. The different lines correspond 
to w = 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, from up to down. 



Figure 2. Convergence of the heavy quark limit S'-wave sum rule (|36|) in the GI 
model [0 for different values of the scaling variable w as one increases N max , the 



number of radial levels included in the truncated variational harmonic oscillator 

N f {n) 

basis. The Ox axis represents (N max + 1) and the Oy-axis j^ r i/2( u 0> that 

— > oo. The different lines correspond 
to w = 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, from up to down. 



the sum rule predicts to be equal to - for iV„ 



19 



IVlUiJri/ij 


M Q 




A/f(0) A/f 

M y '-IVlQ 


m 1/2 -m q 


V M J 


V M h/2 




(GeV) 




(GeV) 


(GeV) 


(GeV 6/2 ) 




ISGW [11] 


10 4 


10 


0.0438 


0.5467 


0.422 


0.235 




10 4 


15 


0.0436 


0.5467 


0.428 


0.235 




10 4 


20 


0.0435 


0.5467 


0.431 


0.236 






infinite 






finite 


finite 


VD [12] 


10 4 


10 


0.119 


0.620 


1.36 


0.603 




10 4 


15 


0.108 


0.620 


1.58 


0.617 




10 4 


20 


0.108 


0.620 


1.76 


0.631 






infinite 






infinite 


finite 


GI [13] 


10 4 


10 


0.386 


0.792 


0.649 


0.620 




io 4 


15 


0.386 


0.792 


0.662 


0.632 




10 


on 

20 


0.386 


0.792 


0.667 


0.640 




10 4 


infinite 






finite 


finite 


CNP [14] 


10 4 


10 


0.389 


0.859 


0.747 


0.669 




10 4 


15 


0.387 


0.858 


0.798 


0.691 




10 4 


20 


0.386 


0.858 


0.828 


0.704 






infinite 






finite 


finite 



Table 1 



radial excitation 


>/M/ (n) 
(GeV 3 / 2 ) 


(GeV 3 / 2 ) 


n = 


0.67(2) 


0.64(2) 


n = 1 


0.73(4) 


0.66(4) 


n = 2 


0.76(5) 


0.71(5) 


n = 3 


0.78(9) 


0.73(8) 


n = 4 


0.80(10) 


0.76(11) 


n = 5 


0.81(17) 


0.77(17) 


n = 6 


0.82(15) 


0.78(15) 


n = 7 


0.82(28) 


0.78(27) 


n = 8 


0.83(25) 


0.79(25) 


n = 9 


0.80(40) 


0.76(40) 


n = 10 


0.83(42) 


0.79(40) 



Table 2 



20 



. 975 



. 95 



0.925 



0.9 



. 875 



0.85 




Figure 1: 



0.5 



0.48 - 



0.46 



0.44- 



. 42 



0.38 



0.36 




Figure 2: 



21 



